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cylindrical surface which is coaxial with the symmetry axis of conical space. The 
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In the case of electromagnetic field, the perfectly conducting boundary conditions or 
isorefractive matching conditions are imposed on the cylindrical surface. For a massless 
scalar field, the semi-transparent conditions (<5-potential) on the cylindrical shell are 
investigated. As a result, the total Casimir energy of electromagnetic field and scalar 
field, per a unit length along the symmetry axis, proves to be finite unlike the case of 
an infinitely thin cosmic string. In these studies the spectral zeta functions are widely 
used. It is shown briefly how to apply this technique for obtaining the asymptotics of 
the relevant thermodynamical functions in the high temperature limit. 
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1. Introduction 

Quantum matter fields on manifolds with conic singularity are of interest in black hole 
physics and especially in cosmic string theory. The status of cosmic strings was changing 
during their rather long history. Therefor it is worth mentioning a few words about this 
topic. 

In the course of time evolution and extension, the Universe should undergo a 
number of phase transitions. From quite general point of view, these transitions may be 
accompanied by emergence of topological defects^ such as monopoles, cosmic strings, 
and domain walls pp. When neglecting the internal structure of these objects, we are left 
with point-like, one- and two-dimensional defects, respectively. For the cosmic scenarios 
the most interesting were monopoles and cosmic strings. On this subject, there is a vast 
literature (see, for instance, [2] and references therein). In the 1980's, the cosmic strings 
were considered as candidates for a mechanism of galaxy formation. However, recent 
observations on the cosmic microwave background have shown conclusively that this 
effect can at most account for a small fraction of the total power, up to 10 % [3J H]. 
Nevertheless, the cosmic strings are still considered as the plausible sources of detectable 
gravitation waves (see, for example, j5j E] and references therein), gamma ray burst [7], 
high-energy cosmic rays [8], and due to their gravitational lensing effects Recently, 
cosmic strings attract a renewed interest in the framework of brane inflation [10] and 
as the macroscopic fundamental strings [TH [121 1131 EH]- Comprehensive review of the 
cosmic string studies are presented in papers [151 EE]- 

The gravitational field, i.e. geometry of the space-time, connected with cosmic 
string, is determined by the solution of the Einstein equations with the energy- 
momentum tensor defined by the string. When infinitely thin straight string with the 
linear mass density /i is located along the z axis then the space outside the string is 
isomorphic to the manifold M 1 x Cm-fa where 1R 1 is an infinite line along the z axis, 
and instead of a plane perpendicular to this axis, one has two dimensional cone C^-^ 
with the angle deficiency <fi = 8nfiG, G being the gravitational constant. The cone 
Cp of an angle becomes a plane when /3 = 2tc. The static straight thin cosmic 
string is thus represented by a locally flat space-time with a conical singularity along 
an axis [TFJ, [TBI EE! ED] . A string with a deficit angle has a positive mass density and it 
is stretched; a string with an excess angle has negative mass density and is squeezed. 
Beside the empty space-time with a single string, cosmic strings appear in a wide variety 
of solutions to Einstein equations. Any axially symmetric space-time can be trivially 
modified to contain a string on the axis of symmetry. Keeping in mind the physical 
origin of a cosmic string, one should characterize it by the gauge flux parameter in 
addition to the deficit angle parameter. However this point is beyond the scope of our 
consideration. 

For physical applications it is interesting to study different fields in the background 

| The topological defects are the regions with higher dynamical symmetry surrounded by a space with 
lower symmetry group. 
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of cosmic string. As explained above, this problem is reduced to consideration of the 
fields on the space-time manifold T 1 x I 1 x C^-^. It is assumed that the effect of the 
fields on the space-time geometry can be neglected. 

One of the difficulties in these studies is nonintegrable singularity at the string axes 
for the energy density of fields considered in the background of infinitely thin cosmic 
string [21] . As a result, the total vacuum energy (per unit length of the string), is 
not defined. To overcome this drawback, the cosmic string of finite thickness can be 
considered. For this purpose one may use, for example, the Gott-Hiscock metric [T8|l22] 
which is the solution to the Einstein equations when the energy density is a constant 
inside the string of a finite thickness and zero outside of it [23] . 

In the present paper, we are going to draw the attention to the following interesting 
fact that is closely related to the problem stated above. Namely, in order to render the 
total field energy in conical space to be finite it is sufficient to impose some boundary 
(or matching) conditions for quantized fields on a cylindrical surface which is coaxial 
with the symmetry axes of the conical space. In a loose form, this observation can be 
interpreted as follows: isolation of a 'naked' conic singularity by introducing additional 
boundary conditions on a surface surrounding it results in 'improving' the physical 
consequences generated by this singularity. 

When choosing the fields for calculating the vacuum energy, we take into 
consideration the following. The vacuum energy of massless fields rule, 
much greater than that of massive fields. Therefore we study the Casimir energy of 
electromagnetic and massless scalar field in the problem in question. As regards the 
boundary or matching conditions on the cylindrical surface, there is a large freedom here. 
We impose such conditions which are usually considered in the analogous studies [241 EH] . 
In the case of electromagnetic field, the perfectly conducting boundary conditions or 
isorefractive matching conditions are widely used in the Casimir calculations. In our 
consideration, we follow in this way. For massless scalar field, we impose the semi- 
transparent conditions on the cylindrical shell. The point is that real boundaries cannot 
confine the physical fields exactly within a given region, i.e. real boundaries are always 
more or less transparent for fields. This property is simulated, for example, by the 
5-potential |25j. 

As a result, the total Casimir energy of the electromagnetic field and the scalar field, 
per a unit length along the symmetry axis, in the conical space with additional boundary 
conditions proves to be finite unlike the case of an infinitely thin cosmic string. In our 
studies of vacuum energy, we use substantially the spectral zeta function technique. In 
addition, this approach enables us to calculate directly the asymptotics of the relevant 
thermodynamical functions in the high temperature limit. 

The layout of the paper is as follows. In section 2, the spectra of an electromagnetic 
field and a scalar massless field in the geometry under consideration are analysed. For 
electromagnetic field, two types of the boundary conditions are considered: perfectly 

§ The manifold M 2 x Cp arises under consideration of fields on the Rindler space-time at finite 
temperature 
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conducting conditions and isorefractive ones. For a massless scalar field, semi- 
transparent matching conditions are applied, i.e. the (^-potential at the point r = a is 
introduced, both signs of the respective coupling constant being investigated in detail. 
In section 3, the global zeta functions for these spectra are constructed by implementing 
the analytical continuation needed. Proceeding from these results, the finite and unique 
values for the vacuum energy of the electromagnetic field and scalar massless field in the 
conical space with additional boundary conditions are derived. In section 4, the high 
temperature asymptotics of the relevant thermodynamical functions in the problem 
under study are obtained by making use of the heat kernel coefficients calculated in our 
previous work [26] and in the present paper. In section 5, we discuss briefly the obtained 
results and compare them with analogous studies. 

2. The excitation spectrum of the fields in conical a space with additional 
boundary conditions 

When the energy-momentum tensor is determined by the 5-like mass distribution along 
the z-axis with the linear density \x the solution to the Einstein equations is given by 
the following metric [HI EB1 H21 12D| : 

ds 2 = -dt 2 + dr 2 + (l-4/iG)Vd^ 2 + d;z 2 , < tp < 2vr , (2.1) 

where G is the gravitational constant. The cylindrical coordinates (r, <p, z) are used 
with the z-axis coinciding with thin straight infinite string. It is easily convinced that 
the space part of metric ( 12. lj) describes the geometry on the surface of a cone with the 
deficit angle $ = or with the opening angle 2n — $. 

When considering matter fields, on the cone surface one can obviously substitute 
the latter by a plane with a segment of the angle $ cut out, provided the periodicity 
conditions in a new angular variable 8 

/(0) = /(0 + 27T-$), O<0<2tt-$ (2.2) 

are imposed. 

In terms of the angular variable 9 metric (12. ip looks exactly like the metric of the 
Minkowski spacetime 

ds 2 = -dt 2 + dr 2 + r 2 d6 2 + dz 2 (2.3) 

supplemented by conditions (12. 2 p (see for details [26J). 

2.1. Electromagnetic field 

In the case of an electromagnetic field in a conical space, we apply perfectly conducting 
conditions or isorefractive matching conditions on an auxiliary cylindrical surface. 
As a result, the problem under consideration is reduced to the calculation of the 
electromagnetic vacuum energy for a cylindrical shell by imposing on its lateral surface 
boundary conditions (or matching conditions) mentioned above and by substituting the 
usual 27r-periodicity by the periodicity condition ( 12. 2D . 
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The solutions to the Maxwell equations with boundary conditions on the circular 
cylindrical shell have been considered in many papers (see, for example, [2T| |26| 127]). 
We shall take advantage of this analysis and allow for the periodicity condition (12. 2(1 . 
The latter requirement leads to the substitution of the integer index n in the Bessel and 
Hankel functions by np, where 
1 $ 

- = 1 . (2.4) 

p 2tt V ; 

According to the physics of cosmic string formation [21 [Til IS] , p is very close to one 
and p > 1 . 

Let us address now the spectrum of electromagnetic oscillations in the configuration 
under consideration. In view of the axial symmetry of the conical space the eigenvalues 
of the spectral problem (eigenfrequencies) u q are 'numbered' by the following three 
indexes: 



which correspond to the coordinates {r, 8,z}. The radial index n r can take discrete 
values or it can be continuous. The angle index n corresponds to the compact angle 
variable 6: < 9 < 2n — $ = 2ii/p. The continuous variable k z in f 1 2 . 5 (1 is the wave 
vector of the plane waves propagating in z-direction: — oo < k z < oo. 

General solution to Maxwell equations in metric (12. 3 p in unbounded space-time 
is expressed in terms of two independent scalar functions that correspond to the TE- 
polarization and TM-polarization of the electromagnetic field [281 EHl ED] ■ In the conical 
space, these functions should satisfy the periodicity conditions (12. 2p . 

For further specification of the spectrum in the problem at hand, the boundary 
conditions on the cylindrical shell should be taken into account explicitly. 

2.2. Perfectly conducting boundary conditions 

On the surface of a perfect conductor, the tangential component of electric field (Ey) 
and the normal component of the magnetic field (Hi) should vanish 



It is important to note that these conditions do not couple the polarizations of 
electromagnetic fields, as the Maxwell equations do, and furthermore the fields inside 
and outside of the shell are independent [3T] . 

Inside the shell, the (unnormalized) eigenfunctions are [26] 




q = {n r ,n, k z } 



(2.5) 



E,| = 0, Hl = 0. 



(2.6) 





n = 0, 1, 2 



(2.7) 



where the parameter p was introduced in (12.4(1 and 
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For simplicity we drop here the common multiplier exp (— iut + ik z z). Solutions to the 
Maxwell equations are expressed in terms of the eigenf unctions ( I2.7P through the well 
known formulas [27J EE] • 

For the TE-polarization, \ nrir in (12. 7p is the n r -th root of the equations 

■4(W*) = 0, n = 0,l,2,..., TV = 1,2,..., (2.9) 

and for the TM-polarization, we should take the roots of the equations 

J np {Kn r o) = 0, n = 0,l,2,..., n r = l,2, .... (2.10) 

The explicit form of the eigenf unctions ( 12.71) implies, specifically, that at n ^ all 
the eigenvalues in the spectral problem under consideration are twice degenerate (the 
both functions svnnpO and cos np8 work) and at n = this degeneracy is eliminated 
because sin np6 vanishes. 

Summarizing we can infer the following: At fixed value of k z the spectrum of the 
frequencies u(TE,k z ) and cu(TM, k z ) inside the shell is discrete and it is given by the 
formula 



cVP + fcf, (2.11) 



where 



w = w nrv (a, k z ), a = TE, TM 



and for the TE-modes, A in (12. lip should be substituted be the roots of the equation 
(12 .9p and in the case of the TM-modes the roots of the equation (I2.10p should be used. 
Outside of the shell the spectrum of the electromagnetic oscillations is continuous 

ck z < lu < oo , (2.12) 

because we are dealing here with free unbounded configuration space. For the 
summation over such a spectrum to be accomplished the scattering formalism should 
be employed [32], i.e. one has to calculate the pertinent S-matrix. The regular solutions 
defining the Jost functions a„ (A) in our model read 

f n (r) = a-(\)H- p (\r) - a+(X)H+(\r), r>a, (2.13) 

where A is a continuous variable related to u and k z by ( 12. 8ft and H^(Xr) are the Hankel 
functions 

H+{Xr) = Hg-\\r), H~ (Ar) = (Ar) . 

The regular solutions ( 12.131) should satisfy the boundary conditions 
d 



alr Jn 1 



= 0, /n™(r)L a = 0. (2.14) 

=a 

By making use of the definition of the S'-matrix 

S n (u) = ^ 

a nW 
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and fl2TT3j) . fl2TTIj) we obtain 

5 " M "^M' 5n (2 ' 15) 

The note concerning the different degeneracy of the states with n = and n / 
for oscillations inside the shell applies also to the scattering states outside of the shell. 

When u 2 < c 2 k 2 there are solutions to the wave equation which decay in the radial 
direction 

2 



-\/ c 2 k\ — to 2 



However they do not meet the boundary conditions (12.141) . Hence in the problem at 
hand there are no surface (or evanescent) waves. 

2.3. Isorefractive (Diaphanous) matching conditions 

When employing these conditions we assume that the electric and magnetic properties 
of the media inside (ei, fii) the shell and outside of it (£2,^2) are different, but the 
velocities of light in these regions are the same: c\ = C2, q = cj y/EiH^ % = 1,2 (c is 
the light velocity in vacuum). The fields E and H inside the shell and outside of it are 
coupled due to the continuity requirement for their tangential components 

discont (E||) = 0, discont (Hy) = 0, r = a. (2.16) 

The configuration space is noncompact, and as a consequence the electromagnetic 
spectrum is continuous (I2.12p . Indeed, for a given n and k z the scattering solutions to the 
Maxwell equations contain 6 amplitude: the TE-solution has 1 internal amplitude and 
2 amplitudes outside of the shell. The same holds for the TM-solution. The matching 
conditions ( 12.161) at the cylinder surface lead to 4 linear homogeneous equations for 
these amplitudes. Hence no restrictions arise here for the spectral parameter u 2 /c 2 . The 
polarizations in this problem decouple in the equations determining the eigenfrequencies 
[33] . but they do not decouple on the level of the scattering solutions!]]] However, the 
scattering matrix in this case has a simple structure (more precisely, the Jost matrices 
have zero diagonal elements). In our paper [32] the S- matrix has been derived for 
the scattering of electromagnetic waves on an infinite circular material cylinder in the 
general case when c\ 7^ c^- These formulae are considerably simplified if c\ = Ci- As 
before, the S- matrix is the (2x2) matrix acting on polarizations (TE and TM) and it 
satisfies the matrix equation^ 

K + S + K- = 0, (2.17) 

|| This is also true for the natural modes of a compact cylinder with c\ = c-i (see (6) and (7) in section 
9.15 of the book [31]). This set of equations does not split into equations separately for a^, and for 
b„, b^, when c% = c 2 . 

It is convenient to interchange the notations for the matrices in [32]: — > K T and take in 
the right-hand side of (20) in [32] the difference instead the sum in order to comply with (|2.13|) in the 
present paper. 
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where the matrices K ± are now 



K 



± 



± 





7 ± 




(2.18) 



rip 



— 1- 



7, 



cJ, 
uX 



up 



cJ, 



up 



— (^2-JnpH^p — /ii J np 'H^\ 
ip > ' 

{^2JnpH np — SiJ np H np S j . 



In these formulae, all the Bessel (J np ) and Hankel (H^ ) functions have the same 
argument aX, where A = a/ (^/c) 2 — k 2 z . Ultimately, the determinant of the S- matrix 

detS = -p4- (2.19) 
det K+ v ' 

is expressed in terms of the known multipliers A™ and A™ [27|,[33l|3l], these multipliers 

being constructed for the both Hankel functions of the first (Hnp) and second (Hnp) 

kinds 

It is easy to show that in the problem under consideration there are no 
surface modes. For u 2 < c 2 k 2 z there are solutions to the wave equation which 
decay in both directions from the lateral surface of the cylindrical shell: u n (r) ~ 
In P {r\Jc 2 kl —u 2 /c), r < a and u n (r) ~ K nv [r ^/ c 2 k 2 z — tu 2 /c), r > a. The respective 
frequency equations are 



I ' 

-'rap 



L np 



£2 



K ' 

1 v np 
- Iv np 



I'l 



1 np 



L np 



fJ>2 



K ' 
Jv np 







n 



0,1,2,... . 



(2.20) 



Here all the modified Bessel functions have the same argument (a/ c)^c 2 k 2 — cu 2 . Since 
lnp > and K np ' < the left-hand sides of Eqs. (I2.20p are strictly positive. Hence 
these equations have no real roots in the interval < u < ck z . 

Summarizing this subsection we infer that the spectrum of electromagnetic 
oscillations in the case of isorefractive matching conditions is pure continuous (see 

(TO). 



2.4- Semitransparent matching conditions for scalar field 

Now we address the consideration of a massless scalar field <p(t, x) in the conical space 
with a cylindrical surface of radius a. On this surface we impose on the field (p(t,x) 
the matching conditions which may be interpreted as the scalar (^-potential. The wave 
equations in this case read 

i o 2 



c 2 dt< 



- A + -5(r - a) 
r 



(p(t, x) = . 



2.21 



In [34] (see endnote 18 there), it was proposed to do this 'by hand' without rigorous justification. In 
■3] rather complicated contours were used when gong to the imaginary frequencies. 
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where g is a dimensionless constant specifying the strength of the potential and A is 
the Laplace operator for the spatal part of metric (I2.3p . Separating the variables 

/ sin npQ s 
\cos np6 y 

we arrive at the equation for the radial function f n (r) 



-iivt+ik z z 



fn(r) 



n 



0, 1, 2, 



(2.22) 



d 2 f n (r) , ldf n (r 



+ 



+ 



2 2 2 

u z 2 n z p z g 

— -K - 2 $( r 



air 2 ' r air 

First we demand the function f n (r) to be continuous at the point r 
fn(a + 0) - f n (a - 0) = 0, 



/ n (r) = 0. (2.23) 



(2.24) 



so that one could integrate the left-hand side of (I2.23P over rdr. Such an integration in 
the viicnity of the point r = a yealds 



9 



/> + 0)-£(a-0) = */„(o). 



(2.25) 



As usually [35] the wave equation with the ^-potential f !2.23|) is substituted by a 
pertinient free equation (without potential) 



al 2 f n (r) , l# n (r) , /w 



<ir 2 



+ 



r dr 



+ 



fc - 



n 2 p 2 



fn(r) = 



(2.26) 



the solutions of the latter, f n (r), being subjected to the matching conditions (I2.24p and 
(I2.25p . Sometimes the matching conditions generated by the (^-potential are refered to 
as the semitransparent ones. 

Finally at a given k z , — oo < k z < oo we face the spectral problem specified by 
differential equation (12. 26ft . matching conditions (I2.24p . (12. 25ft and physical conditions 
at the origin (r = 0) and at the infinity (r — > oo), the frequency u > being the 
spectral parameter. It turns out that the structure of the spectrum in the problem at 
hand depends strongly on the sign of the constant g, namely, for positive values of g the 
spectrum is pure continuous, and when g < we have, in addition to the continuous 
branch, surface modes (or bound states). 

First we consider the continuous part of the spectrum, i.e. we construct the 
scattering matrix in this problem for both signs of g. The regular solution is 



fn(r) 



a™(\)J np (\r) 



for r < a, 
for r > a , 



where 



A = v^/c) 2 - kl . 
The matching conditions (12.241) and (I2.25P yield 

a™(\)Jn P (\a) = a~{\)H~ p {\a) - a+{\)H+ p {\a) , 

C(A) [j-JnAXa] + J np '(Xa)] = a-(X)H- p '(Xa) - a+(\)H+ p '(\a 



(2.27) 

(2.28) 

(2.29) 
(2.30) 
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Here and below the prime in Bessel and Hankel functions denotes the differentiation 
with respect to their arguments, Xa. Eliminating from f |2.29j) and ( 12 .30 j) the amplitudes 
a^ n (A), we obtain the S'-matrix 

gJn P (Xa)H- p (\a) + 2i/n 
n[U) ~ gJ np (\a)H+ p (\a) - 2I/7T- ^ 

Now we address the consideration of the surface modes in the problem under study. 
Such solutions to the boundary- value problem (12.261) . ( 12.2411 . (12.251) exists only for 
negative constant g. Indeed, when 

— ~k 2 z <0, 
cr 

the radial wave equation 02.26(1 has the solutions K np (K,r) and J np (/tr), which 
exponentially vanishe outside the cylinder shell in both directions, for r > a and for 
r < a respectively. Here 



K=^kl-^. (2.32) 

Thus, the complete solution describing the evanescent waves (surface modes) is given 
by 

{a™(K,)I n JK,r) for r < a, 

' (2.33) 
a^(K)K np (nr) for r > a. 

The matching conditions (12.241) and (I2.25P give 

<»/ np M = a e n x (n)K np (Ka) , (2.34) 

a™{K)K' np {Ka) - a?:{K)I' np {Ka) = ^-a^I^a) . (2.35) 

CLKj 

Taking into account the relation [36] 

I„{x)K' v {x) - I' u (x)K„(x) = -- , 

X 

we obtain equation determining the frequencies of the surface modes: 

1 + gI np (x)K np {x) = , (2.36) 

where 



x 2 = aV = a 2 (k 2 z - ^ j . (2.37) 

The product of modified Bessel functions I v {x)K v (x) has the following properties 
I (x)K (x) ~ lnx, x — > , 

7,(0)^(0)=^, u>0, (2.38) 

Au 2 - 1 



I v {x)K v (x) ~ ^- 
2x 



(2 X y 



+ 



x — > oo . 
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(see, for example, [US])- By numerical calculation^ one can easily make sure that, for all 
u, the product I u (x)K u (x) monotonically decreases in the interval < x < oo varying 
between the limiting values given in (I2.38p . Hence, the frequency equation (12.361) with 
g < has always a root for n = but for n > 1 this equation has the root, provided 
that the inequality 

W"* < £ (2-39) 

is satisfied. Obviously for a given \g\ the condition (I2.39|) can be true only for some first 
values of n. As a consequence, merely for such n the frequency equation (12.361) can have 
solution. 

Let x — clk be a root of f l 2 . 3 6 j) . then the frequency of the surface mode, oo sm , is 
given by 

Wsm = C^kl - K 2 , (2.40) 

(see (I2.37P ). Obviously, the frequency w sm should be a real quantity thus one has to add 
to (I2.40p the restriction 

k z >n. (2.41) 

Summarizing we can claim that the spectral problem under consideration has the 
continuous spectrum in the domain f 1 2 . 1 2 j) for both signs of the coupling constant g. For 
negative g there appear, in addition, the surface modes with frequencies (12. 40 p . (I2.4ip . 

Closing this section it is worth mentioning solutions in the problem under 
consideration which are obtained from the surface modes when k z < k. These solutions, 
as the surface modes, are located around r = a (see (I2.33P ). but they have pure imaginary 
frequencies 

u) = ±iu with co = cy k 2 — k 2 , k > k z . (2.42) 

As a result, such solutions exponentially rise with time ~ exp(±ut). Hence, already at 
the first-quantized level, there is instability in the system at hand. However, in what 
follows we disregard this instability. In our paper [37] analogous solution was discussed 
briefly in the framework of the flat plasma sheet model with negative parameter q, 
and it was interpreted as a resonance solution. In reality, it is not the case and it is 
impossible to confront the considered solutions of the Klein-Gordon equation with any 
special solution of the Schrodinger equation. 

3. Global zeta functions and Casimir energies 

Now we are in position to construct the spectral zeta functions 

c(*) = EK 2s -^~ 2s ) 



* The authors are indebted to M. Bordag for such calculations. 
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and to calculate, on this basis, the vacuum energy of quantum fields in question 

E=U(s = -\). (3.2) 



2 V 2 , 

The summation in (13.11) should be done over the whole spectrum (discrete, continuous, 
and with allowance for the surface modes if such exist) . The frequencies u q in (13.11) are 
obtained from u q when a — > oo. The subtraction of uj q in (13.1 p corresponds to removing 
from the vacuum energy the contributions proportional to the volume of the manifold 
and to the manifold boundary In this way the renormalization of the vacuum energy 
is accomplished. The parameter s is considered at first to belong to the region of the 
complex plane s where the sum in (13.11) exists. After that the analytical continuation of 
(13. ip to the point s = — 1/2 and further for Re s < —1/2 should be done. When summing 
over the continuous branch of the spectrum, the function of the spectral density shift 
have to be employed [32]. The rigorous mathematical theory of scattering gives the 
following expression for this function 

Ap(u) = p(u) - po(u) = t^^j det S(u) , (3.3) 

where S(u) is the S"-matrix in the spectral problem at hand. Here p(k) is the density of 
states for a given potential (or for a given boundary conditions in the case of compound 
media) and po(k) is the spectral density in the respective free spectral problem (for 
vanishing potential or for homogeneous unbounded space). 

3.1. Perfectly conducting boundary conditions 



As shown in section |2T2| the spectrum in this problem has two branches: discrete branch 
(electromagnetic oscillations inside the shell) and continuous branch (oscillations outside 
of the shell) . In a complete form (13.ip reads now 



00 oo 

C( S ) = 2 J2 / ^E' Y,"nnM,k z )+ [u-^Ap n (a^,k z )du 

a 71=0 n r , 

— oo L ck z 

where a = TE, TM and the prime over the sum sign means that the term with n = 
is taken with the weight |. The discrete frequencies u nnr (a, k z ) in (I3.4p are defined by 
( 12.91) and ( 12.101) . and the function of the spectral density shift is given by ( 13. 31) and 

The further treatment of (13 .4p is aimed at deriving a unique integral representation 
for both contributions of discrete branch and continuous branch of the spectrum, the 
integration being done over imaginary frequencies u) — > iu. This is accomplished by 
transforming the sum Y2 n i^o the contour integral applying the argument principle 
theorem from the complex analysis (see figure la). After that the initial contour C is 
continuously transformed into the contour Cq with Q — y oo. The contribution of the 
continuous part of the spectrum is also expressed in terms of integrals along the contours 
Cq shown in figure (lb). The contour (Cq) is used for transforming the integral 
f°£ dujuj~ 2s ... in (13.41) containing H + (H~) and H + ' (H~') in accord with (12. 15ft . In 



(a -> oo), (3.4) 
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Figure 1. The contours on the complex oj plane which are used when going on to the 
integration over the imaginary frequencies in the case of bounded (a) and unbounded 
(b) configuration spaces. 



both the cases the integration is accomplished on the complex frequency plane cu with 
the cut connecting the points ck z and — ck z . 

It is worth noting here that the analytic properties of the scattering matrix (or 
the Jost function) in the Casimir studies are different in comparison with the standard 
theory of potential scattering |38J. In fact they are close to those for the Klein-Gordon 
equation, the role of mass squared in Casimir calculations being played by k 2 z . This 
implies in particular that the analytical properties of the scattering matrix in the Casimir 
calculations should be revealed by a direct analysis of its explicit form without referring 
to the nonrelativistic potential scattering. 

Taking into account all this we obtain the following formula for the spectral zeta 
function in the problem under consideration: 

OO 

v/?ar(s)r f - - .* i 



where 

Vn P (y) = y [in P (y)K np (y)]' . (3.6) 

It is worth comparing (13. 5p with (3.11) in [31]. Further we will follow the well elaborated 
procedure for construction of the analytical continuation of the formulae like (13.51) into 
the left semi-plane of the complex variable s (see, for example, [39J). This regular 
procedure substantially uses the uniform asymptotic expansion (UAE) of the modified 
Bessel functions. Remarkably, the result obtained in this way can be of arbitrary 
required accuracy. To obtain it, one should keep in UAE sufficiently many terms. As 
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far as we know, there are no other methods to construct the analytical continuation 
needed. 

Upon the change of the integration variable in (13.51) y = npz, n = 1, 2 ... we can 
use here UAE. We content ourselves with the first two terms in UAE 

ln<!l 



zH* 



An 2 p 2 



1 + 



An 2 p 2 



z—(I np (npz)K np (npz)) 



3 - 30t 2 + 35t 4 + ^ z 4 t 4 + 0(n" 4 ) 



2N-1/2 



(3.7) 



Now we substitute (13.71) into all the terms in (13. 5 \} with n ^ 0. The term with n = 
in this sum will be treated by subtracting and adding to the logarithmic function the 
quantity 

1 y 4 



4(1 + ^)3- ( 3 ' 8 ' 
As a result, the zeta function Csheii(s) can be presented now as the sum of three terms 



Cshell(s) = Zi(s) + Z 2 {s) + Z 3 (s) 



where 

ZAs) - 



2s - 1 



Zo(s) 



2^aT(s)T [~-s 



c\~2s (l-2s)(3-2s) 



-2s 



dyy 



-2s 



Hi -ti{y)] + - A y A t\y) 



(3.9) 
(3.10) 



r 



[2 P - 1 ^ s Cr(2s + 1) + 1]- 



r( s ) 



c\ -2* (1 - 2s)(3 - 2s)(284s 2 - 104s - 235) 



r [- + S 



r s) 



(3.11) 



Cr(3 + 2s). (3.12) 



6 1440 y^ap 3 + 2s 

As was shown in [31], Zi(s) is an analytic function of the complex variable s in the region 
—3/2 < Re s < 1/2. In (13. lip and (13.121) Cr( s ) is the Riemann zeta function which 
accomplishes the analytical continuation of the sum Y^=i n ~ 2s t° the left semi-plane 
of s. 

It is left now to take the limit s — > —1/2 in (13. 9p . Contribution Zi(s) (I3.10p to the 
zeta function does not contain the parameter p, i.e. it remains the same as in the case 
of Minkowski space-time. Hence for Z\{— 1/2) the value from [31] can be taken 



Z 1 (-l/2) 



2vra 2 



;-0.6517) 



(3.13) 



A special care should be paid when calculating the limit s — > —1/2 in Z 2 (s) (13.111) in 
view of the poles of the function T(s + 1/2) at this point. Using the values 



Cr(0) = -~ C R (0) = -^ln27r, T(x) = 1 - 7 + 0(x) 



(3.14) 
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we derive 



lim [2 p - 1 - 2 Xr(1 + 2s) + 1}t(1 + s 

s— >— 1/2 \Z 



= lim {2[1 - (1 + 2a) hip + . . .][Cr(0) + Cr(0)(1 + 2s) + ...] + 1} - 7 + . . 

8-t—l/2 \L + ZS J 

2tt 

= -2 In— . (3.15) 
P 

Thus we have for Z 2 (— 1/2) 

Z 2 (-l/2) = -^ln-. (3.16) 
o7rer p 

The third term (I3.12p gives 

Z 3 (-l/2) = -At — ^ • (3.17) 
3V 7 ; 2vra 2 480 6p 2 V ; 

Gathering together Zj(— 1/2), i = 1, 2, 3 we obtain the finite value for the vacuum 
energy of electromagnetic field in the problem under consideration: 

£.„.,, = !c. M1 (-i), (3.18) 

where 

/ 1\ C ( 1 27T 7 7T 2 1 \ . . 

("2 j = w I" ' 6517 + 1 ln 7 + «oT?J ' (319) 

At p = 1 we reproduce (3.35) from [53] . 

The consideration presented here can be extended to the next terms in the UAE of 
the Bessel functions (13 .TP in a straightforward way. Therefore we shall not present here 
these rather cumbersome expressions. 

3.2. Diaphanous matching conditions 

In this case, we have continuous spectrum of electromagnetic oscillations (I2.12|) . As 
a result only the second term in (13. 4p survives. When going to the imaginary 
frequencies the contours should be used (see figure lb). We again obtain the integral 
representation (13 .5p for the spectral zeta function, but now fi np (y) should be multiplied 
by the parameter £ 2 : 

e = = * I- (3-20) 

Performing in the same way as for the perfectly conducting boundary conditions, we 

arrive at the result 

*n„ M\ < 2 / 1 2tt tt 2 tt 2 7 1 \ , N 

d 1 " — = _A_^ -0.490878 + - ln — + . 3.21 

Scyl V 2 J 2na 2 \ 4 p 288 360 1920 p 4 / 

We are considering only the contribution to the zeta function which is linear in £ 2 . 
When p — y 1 the derived value for (^(—1/2) vanishes [331 SU S2]. As was noted earlier, 
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in cosmic string physics the parameter p is close to 1. Hence, the vacuum energy of 
electromagnetic field for these boundary conditions 




(3.22) 



will be also small. 



3.3. Semitransparent matching conditions 

At first, we assume that g > and, as a consequence, the spectrum of excitations of the 
field <p(t,x) is pure continuous (see section l2Ui) . When constructing the spectral zeta 
function, we carry out the integration over the continuous spectrum (I2.12p by making 
use of (13. 4ft keeping there the second integral term alone. Further there is no need to 
subtract here the contribution of Minkowski space-time (the term obtained in the limit 
a — > oo). The point is that such a subtraction has been already done in ( 13. 31) . By 
making use of the exact expression for the ^-matrix ( 12.311) and the integration contours 
shown in figure (lb), we obtain the following representation for the zeta function 

Cs-t(s) = C(s) ]T' dyy^^il + glMKM] , (3-23) 

n=0 J ° ^ 



where v = np and 



c -2s 



Further we shall consider the zeta function f l3.23f) in the linear approximation with 
respect to g 

( s - t (s)=gC(s)f^ r dyy^^[lMKM] • (3-25) 

n=0 ^° ^ 

Analytical continuation of this function to the left semi-plane of the complex variable s 
will be accomplished in the standard way (see, for example, [39J): we add and subtract 
to the integrand in (I3.25P for n > 1 a few first terms of the UAE for the product of the 
modified Bessel functions l v {yz)K v {yz): 

I v (vz)K v (vz)*£^, t = ^=, (3.26) 



where 



k=0 

P (t) = | (3.27a) 

^(^ = ^(1-6^ + 5^), (3.276) 
16 

t 5 

P 2 (t) = (27 - 580t 2 + 2170t 4 - 2772t 6 + 1155t 8 ), (3.27c) 

256 

+7 

P 3 (t) = 7^(t 2 - l)(425425t 10 - 1106105t 8 + 1014442t 6 (3.27 d) 

- 383570t 4 + 51445t 2 - 1125) , (3.27e) 
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and so on 

Do not pretending to high accuracy, we take in UAE ( I3.26P for subtracting and 
adding only two terms 



I v {yz)K v {yz) ~ ^ 



P k (t) 



k=0 



2/fc+l 



t 

2^ 



t 2 (l - 6t 2 + 5t 4 ) 

1 + — " + 



(3.28) 



In the case of the product I (y)K (y) we add and subtract the usual asymptotic of this 
expression when y — > oo 
1 



, =. (3.29) 

2v^TV 

As a result, the zeta function (13. 25j) can be represented now as the sum of four 
terms 

4 



i=i 



where 
Z x {s) - 

Z 2 (s) - 

Us) - 



l-2s 



m 

\n=l 
1 oo 

16 ^ 



v 



dy 
2 



Ia(y)K (y) 



2^1 + 2f 



dzz L - 2s —t 
dz 



-2-2s 



n=l 



dz Z 



l-2s 



r(i - 6r + 5f 



oo „ 
n=l J ° 



oo T 
l-2s u 



dz z 



dz 



dz 

I v {vz)K v (vz) 



t 

27 



1 + 



t 2 (l - 6t 2 + 5t 4 



(3.30) 

(3.31) 
(3.32) 
(3.33) 
(3.34) 



The obtained formulae (see (I3.3ip . (I3.32p . and (13.331) ) show clear the basic idea of 
making use of the UAE for analytic continuation: this method allows one to factorize 
the divergences which are originated in the sum over n and in integration over y or z. 
In what follows we express the result of the summation in terms of the Riemannian 
zeta function and the results of integration through the gamma function. It is this that 
provides us with explicit analytic continuation needed. The last term, Z^s), can be 
evaluated only numerically upon fixing the parameter p [y = np) and the value of the 
variable s, for example, s = —1/2. As a rule, this term gives a small correction [33| 140]. 
Fortunately, in the problem under consideration this term seems do not contribute to 
the Casimir energy for arbitrary values of p (see below). 

Taking into account the behaviour of (d/dy)[Io{y)Ko(y)] at the origin and at infinity 

1 

— , y^o, 
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we infer, that the function Z\(s) defined by (I3.3ip is an analytic function in the region 

-KRes<l/2. (3.36) 

The upper bound in ( 13.36!) > which is caused by the behaviour of the integrand in (I3.3ip , 
can be removed by introducing the infrared cutoff or the mass of the field </?(£, x) (see, 
for example, [ffij p. 4524]). At the end of calculations this cutoff (or mass) should be 
set zero. At the point s = — 1/2 we obtain numerically 

g C(-l/2)Z 1 (-l/2) = • 0.999 .... (3.37) 
Integral in ( 13 .32 j) exists originally in the region 

< Re s < 3/2 , (3.38) 
and the sum over n in this formula is finite in the domain 

Re s > 1/2 . (3.39) 
Thus, in the region 

1/2 < Re s < 3/2 (3.40) 

Equation (I3.32p defines an analytic function of the variable s. We can analytically 
continue this function to the whole complex plane s, except for simple poles at some 
isolated points. For this purpose we have to put 

oo 

£ v~ 2s = P- 2 Xr(2s) , (3.41) 



n=l 



-i-a. d + f°^2-2 S+ 3 i r( a )r(3/2 - a) 



dz z'- 23 -^- 1 = - dz z 2 - 2s t d = v ' . (3.42) 

dz 7 2 T(3/2) 1 ; 



We have used here the table integral 

= i V . (3.43) 

The right-hand side of (13.431) has no singularities (poles) in the region 

Re a < 0, Re (a + p) > . (3.44) 

Admitting the appearance of poles, we define the left-hand side of (I3.42p in the whole 
complex plane s by the right-hand side in this formula. 
Finally we have 

gC(s)Z 2 (s) = - 9 -C(s) [2p-*&(2«) + l] ^ fj*~ $) 



gc 2s 



2na 



l-2s 



(3.45) 
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At the point s = —1/2 it gives 

(l- |). (3.46) 



gC(-l/2)Z 2 (-l/2)-. ''" 



Kr(-i) + ~ 



eg 

Ana 2 



2wa 2 

In the same way we can construct the analytic continuation of the function Z^(s), 
which is defined originally by ( I3.33P in the region 

0<Res<3/2. (3.47) 

In the whole complex plane s, we have 

g C(s)Z 3 (s) = - ^W-^ Cr(2 + 2S){2S + 1)(2S ~ 1} ' (3 ' 48) 
At the point s = — 1/2 we find 

g C(-l/2)Z, { -l/ 2 ) = -£L- p . (3.49) 

The functions Zi(s), i = 1, 2, 3 are defined by initial formulae (I3.3ip - (I3.33f) in the 
domains of complex plane s which have a common strip, namely 

l/2<Res<3/2 (3.50) 

(see (13.36P with the note following, (I3.40p . and (I3.47P ). This implies, that the analytic 
continuations of individual functions Z^s) i = 1, 2, 3, constructed above, applies also to 

the sum J2i=i z i( s )- 

Gathering together the contributions (I3.37p . (I3.46p . and (I3.49j) we get 

C_ t (-l/2) = -^fp-iV (3.51) 



247ra 2 \ p / 

It seems that this answer is exact (in the linear in g approximation). The point 
is the last term, Z 4 (s) in the sum (13 .30 j) does not contribute to the vacuum energy for 
arbitrary p. This assertion is based on the following fact noted, for the first time, in the 
paper [HJ Appendix A]: Upon substituting in (13. 34[) the product l v {yz)K v {yz) by the 
uniform asymptotic expansion (13.261) and setting here s = — 1/2, we obtain 

°° poo J 

= W- 2 *Cr(2A;-1) / dzz 2 -^-P k (t). (3.52) 

k=2 J ° ^ 

Analytical integration in (I3.52p gives zero 

dzz 2 ^-P k (t) = 0. (3.53) 
az 

We have checked this relation for k = 2, 3 by making use of the table integral [13] 

d T ( dL ^ A ) T (p-^ A ) 

dzz l - s -^t 2{p - 1] = (1-p) V / —\ 3-2Rep<Res<3. (3.54) 

dz T{p) 
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Certainly, it is important to bring to light the origin of these null results (see also 
[251 E]). 

Proceeding from all this, we conclude that the exact value (in the linear in g 
approximation) of the Casimir energy in the problem at hand is 

£U = ^(-l/2) = ^( P -i), 9 >0. (3.55) 

The vacuum energy ( 13.551) is positive because we have assumed that the constant g 
is positive. In the limit p = 1 (13.551) reproduces the vanishing vacuum energy of the 
cylindric delta-potential (25j |44] . 




Figure 2. Contours in the complex w plane which should be used when going on to 
the integration over imaginary frequencies in the case of spectrum having continuous 
branch and surface modes. There are two cuts along the real axis: — oo < lu < and 

ck z < u) < oo. 



Now we address the consideration of negative constant g, when the spectrum of the 
field (p(t, x) has a continuous part (" 12 . 1 2[) and surface modes with frequencies (12.40)) . In 
conformity with this we represent the zeta-regularized vacuum energy E(s) as the sum 
of two terms 

E(s) = E sm (s) + E cmt {s) (3.56) 

(see ( 13 Ah ). Transition to the imaginary frequencies in E cont (s) is accomplished by 
making use of the integration along the contours Cq and shown in figure 2. One 
can see from this figure, that E cont (s) involves, with an opposite sign, the contributions 
of the surface modes u sm (1 12. 40 ft with \k z \ > k) and the contribution of pure imaginary 
frequencies ±i£j (see ( 12.42)) ) . The contribution of co sm to E cont (s) is canceled with E sm (s) 
(integration along the contour C £ , the contributions of the imaginary frequencies ±ia) 
are mutually canceled at s = —1/2. 
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As a result, we arrive again at the formula (I3.23P with an alone change, namely, the 
integral over imaginary frequencie^l y should be treated as the principal value integral 

E(s) = ^C(s) f^P.V. f™ &yy x - 2s ^-\n[l-\g\I v {y)K v {y)] , (3-57) 

n=0 ^° ^ 

where C(s) is defined, as before, in (13 .24ft . Obviously, in the linear in — \g\ approximation 
we do not 'feel' the P.V.-prescription. Hence in this approximation we have the result 
(13.551) with the negative sign: 

^-t = -^4fp--\ 9 <0. (3.58) 

487TGT \ p J 

In [25j HI] dealing with a massless scalar field in the background of semitransparent 
cylindrical shell (p = 1), the negative values of g were not considered at all with the aim 
'to avoid the appearance of negative eigenfrequencies' [25l p. 29]. The vacuum energy 
of massive scalar field in the background of a cylindrical ^-potential was considered in 
[4"5| H6] for both repulsive and attractive potentials. The possibility of bound states 
were not investigated there and it was noted that for g < the renormalized Casimir 
energy acquires complex values. 



4. High temperature expansions 

Quantum field theory predicts a specific temperature dependence for the energy of 
radiation connected with a heated body. One cannot anticipate here the Planck 
spectrum and the Boltzmann law. These laws hold for the photons being in unbounded 
space and nevertheless possessing nonzero temperature. For such photons a notion 
of black body radiation has been introduced as far back as in pre-quantum physics 
time. It is clear that in the general case, the spectrum of a heat radiation and its 
temperature behaviour should depend on the shape of a heated body. By the way 
this dependence is important for many radio engineering equipments (thermal noise of 
antennae, waveguides and so on). 

For obtaining the asymptotic behaviour of the thermodynamical functions, a 
powerful method of the zeta function technique and the heat kernel expansion [17] 
can be used. It is important that for these calculations it is sufficient to find the zeta 
function and the heat kernel coefficients not for a complete operator, determining the 
quadratic action of the quantum field under consideration, but only for the space part 
of this operator, i.e. one has to calculate first the zeta function that is used for the 
Casimir calculations at zero temperature. This is an essential merit of this approach. 
It is these spectral zeta functions that have been constructed in the preceding sections. 
Obviously, these results can be used for obtaining the high-temperature expansions for 
the thermodynamic potentials in the problems under consideration. 

jj The variable y is not equal exactly to iw, more precisely, it is related to iX in (|2.8|) . 
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The high-temperature asymptotics for the Helmholtz free energy has the form 
[391 HZ] 

F(T) " - 2 C (0) + 90 " Bl/2 4^ Cr(3) - 24 T + (4^ T 



5 > n 



16tt 2 



B 5/2 h 2 

(4tt) 3 / 2 24T 



R / ft \ 2 ™~ 3 

- T E(^yb(^J r(n-3/2)C R (2n-3) ; T oo. (4.1) 

Here 7 is the Euler constant and -B n /2, w = 0, 1,2, . . . are the heat kernel coefficients 
which can be calculated by making use of the spectral zeta function in the problem at 
hand 

-Bn/2 / n-3" 



(47r)3/2 - s hm i ^ + ^-jr( S )C(,), n = 0,l,2,.... (4.2) 

In order to use this formula the function ((s) should be known in the vicinity of the 
following points 

S = I 1,1,0,.... (4.3) 

The heat kernel coefficient B n / 2 is different from zero only when the product r(s)£(s) 
has a (simple) pole at the point s = (3 — n)/2, n = 0, 1, 2, . . . . 

The asymptotic expansions for the internal energy U(T) and the entropy S(T) are 
deduced from (14. ip through the thermodynamic relations 



U{T)= -T 2 ^{T- l F{T)), (4.4) 

BF 

S(T) = T- 1 (U(T) - F(T)) = -—. (4.5) 

For compact manifolds B ~ V and Bin ~ S, where V is the volume of the 
manifold and S is the area of its boundary. In view of this the second term in the 
right-hand side of (14. ip gives the Stefan-Boltzmann law (U(T) ~ VT A ). 

In the case of noncompact manifolds the relative spectral functions (zeta function 
and heat kernel) are employed (see, for example, (HI H9]). When defining these 
functions the contribution due to free space is subtracted. We have noted earlier that 
the scattering formalism implements this automatically. As a result, in our approach 
B /V -> 0, when V -> 00 and B 1/2 /S -> 0, when S -> 00 (see for details [261 P- 446]). 

Summarizing, we may infer, that the approach applied gives in fact the corrections 
to the Stefan-Boltzmann law which are determined by the specific spectrum of 
electromagnetic oscillations (or other fields), i.e. by the geometry of the problem under 
consideration. It should be noted that these corrections are not proportional to the 
volume V of the system, hence we cannot say about the corrections to the energy density 
of electromagnetic energy. The area of practical employment of these thermodynamical 
asymtotics is calculation of the temperature dependence of respective Casimir forces. 
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The heat kernel coefficients needed in (14.1 j) were calculated in our previous article 
[26J and in Appendix A and Appendix B of the present paper. The value of C'(0) can be 
found by technique employed for this purpose, for example, in jl?]. By making use of all 
this we obtain the thermodynamical asymptotics in the geometry under consideration. 

For electromagnetic field with perfectly conducting boundary conditions these 
asymptotics are 

153 1 h 3 c 2 



F(T) ~ _-C'(0) + — -ln- 
v ; 2 s w Q4p a T 2048 4vr 24Ta 3 ' 

U(T) ~ — T 153 1 * V 



16 Ana 

Sir) ~ ic'(0) 4 



2048 4vr 12Ta 3 
3 



64pT 



153 1 h 2 c 2 



11 T / 2048 4tt 24T 2 a 3 



where 



C'(o) = - 

a 



0.53490 + — f 3 In — + 37 - 4 

32p V 2pc f 



47 



For isorefr active matching conditions we have 

F(T) 



|c«>) 



2 7 " W7r3 + Aii T i n 4 



45 ft 3 



f/(T): 

S(T) : 
where 



T 4 27r 3 a^ 3 



h 3 15 



+ 



16 4pa ~~ T 
45 



64p a 1024 A8pTa 3 ' 



-C'(o) - 

2 s v ; 45 



7T 3 a 2 £ 2 T 3 



r ±e_ 

c 3 h 3 16 4pa 



In 



25Qp 2 

45 £ 2 c 2 ft 2 
1024 pa 3 96T' 

45 e 2 c 2 h 2 

Tj 1024 pa 3 94T 2 



Cr(3) 



£ 2 

C'(o) = - 

a 



0.28428 + — (3 In — + 37 
32p V 2pc 1 



4 



27 



128p : 



:Cr(3) 



And finally for a massless scalar field with (^-potential (g > 0) we obtain 

ngT 2 



where 



F(T) 
U(T) 
S(T) 

C'(o) = 



2 ^ K ' h 3 90c 3 



T 4 7T 3 

ft 3 c 3 30 



2 ^ ^T 2 
^ + 12^ 
T 3 27r 3 ^a 2 



1 



12frc 

r 



1 

1 - - 

p 



2 <,( °> + » 45C 3 



P 

6ftc 



_0_ 

7ra 



-0.9818 + 



1 

t; + - In 

144p 2 2 



i-i 

p 



2ty 
P 



(4.6) 
(4.7) 
(4.8) 

(4.9) 

4.10) 
4.11) 
4.12) 

4.13) 

4.14) 
4.15) 
4.16) 

4.17) 



5. Conclusions 



We have demonstrated that enclosing the conical singularity with a cylindrical surface 
and imposing here appropriate boundary conditions on quantized fields render the total 
vacuum energy of these fields finite, at least, in the zeta renormalization technique. 
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It is worth noting here that usually conical singularity is considered as the 
singularity of the curvature tensor [201 ISO]- We have treated this singularity as a non- 
smoothness of the boundary |26j . 

Calculation of the Casimir energy for configuration close to those considered in 
our paper has attracted much attention in recent studies. Ellingsen, Brevik, and Milton 
[5T| |52| |53| |54"1 [55] calculated the vacuum energy for wedge and cylinder geometry which 
includes, as a special case, our cone configuration. The periodic boundary conditions in 
angular variable should be used for this choice. Only this case has rigorous justification 
for calculations conducted (vanishing of the heat kernel coefficient B 2 ). 

The local energy density for a wedge with a coaxial cylindrical shell was considered 
in papers [561 E3 Ell EH [60]. The singular behaviour of this quantity at r = was noted 
in these studies. This singularity forbade the calculation of the respective total vacuum 
energy. 

As a byproduct we have presented here a complete analysis of the spectral problem 
for a scalar massless field in the background of the ^-potential located at a cylinder 
lateral surface, both the 5-wall and the 5-well being considered. 
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Appendix A. The heat kernel coefficients for isorefractive boundary 
conditions 

Here we calculate the heat kernel coefficients for electromagnetic field obeying to the 
isorefractive matching conditions on the cylindrical surface of radius a in the conical 
space. We shall use the relation (I4.2p between the heat kernel coefficients and the 
relative spectral zeta function ((s). As was noted in section [3721 the spectral function 
in the problem at hand is given by (13. 51) . where (J> np (y) should be multiplied now by the 
parameter £ 2 given in (I3.20p . Thus we have 



where v = np, fi u (y) is explained in ( I3.6p . and C(s) is defined in ( I3.24p . Further we 
separate in (lA.ip the contribution due to the term with n = denoting the rest by C( s ) : 




(A.l) 



In the £ 2 - 




(A.2) 



DO 




(A.3) 
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The integral over y exists in the strip 
- 1/2 < Res < 3/2, 

because 
dy 
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(A.4) 



4(ln^ + 7 + ^ )y+0(y s ), y ->• 0; 



d/4 
dy 



2y 3 Ay 1 



+0(y- 7 ), y^oo.(A.5) 



In view of ( 14 .2D it implies that Co( s ) does not contribute to the heat kernel coefficients 
B k with k = 1/2, 1, 3/2. 

The analytic continuation of Co( s ) given in (1A.3I) to the right of the strip (1A.4D . i.e. 
to the region Res > 3/2, is accomplished by addition and subtraction of the y — > 
asymptote 



{oo 



l-2s 



djPo(y)-4(ln5+ 7 +o I* 



SO 



sing/ 



2e 2 c(s) 



(3 - 2s) 2 



ln2-7 




(A.6) 



(A.7) 



With making use of this result we derive the contribution to the heat kernel coefficient 
B generated by Co"?( s ) : 



(4vr) 3 / 2 res T(s)C a 4(s) = -4< 2 



s=3/2 



(A. 



It should be noted that this contributior lffl is due to infrared singularity of the spectral 
density in the problem at hand. 

The analytic continuation of the integral in ( ]A.3j) to the left of the strip (1A.4|) . i.e. 
to the region Res < —1/2, is accomplished by addition and subtraction of the y — > oo 
asymptote (IA.5I) : 
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e 2 



f dyy l - 2s ^l{y)+ 
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l-2s 



d 2 _ 1 



dy 



+ 2y* ' V 



+ Co7(s) 
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So 



sing/ 
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+ t: 



21 + 2s 4 2s + 3 



(A.9) 
(A.10) 



From ( 1A.10I) we derive the respective contributions to the heat kernel coefficients B2, 
B5/2, B 3 : 



o c 3 o c 

<S' 0, -vr£ 2 - 
cr 4 a 4 



(A.11; 



tt I n EI] the calculations analogous to (|A.6|) and (IA.7I) were done incorrectly. This resulted in erroneous 
conclusion drawn in paper |47j , that Bq = for a material cylinder with c\ = C2 . 
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Now we turn to ((s). To draw out the singularities in ((s) we apply the uniform 
asymptotic expansion of the modified Bessel functions. In the ^-approximation it gives 



ln[l-^H] = -f 



,zH 6 

Au 2 



1 + -r^(3-30t 2 + 35t 4 ) 



(A.12) 



t 4 



+ —(9 - 256t 2 + 1290 1 4 - 2037 1 6 + 1015 1 8 ) 
4z/ 4 ' 



The integrals over z converge in the range (5 —j)/2 < Res < 5/2, where j is the power 
of t in (IA.12I) . After the integration over z in the divergent parts of the zeta function 
£(s) we arrive at the result 

■■2_2s-l„l-2«-i 



-C R (2s-l + z) Ai(s) 



20Fc 2s r(s)r(3/2 - s) 

A 2 (s) = -— , (-1 + 2s) 2 (-3 + 2s), 



1,2,..., 



(A.13) 
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1 7T 



24576 cos(7rs y 

1 7T 



(-3 + 2s)(27 - 46s - 44s 2 + 56s 3 )(-l + 4s 2 ) 
(-9 + 4s 2 )(-l + 4s 2 ) 



47185920 cos(vrs) 
(9586s + 4640s 5 - 18000s 3 + 15272s 2 - 2576s 4 - 6345), 

A 2i+ i = 0, 1 = 1,2,...,. 
For obtaining ( ]A.13j) the integration was carried out with making use of the formula 

dm. 

The contributions of £(s) to the heat kernel coefficients Bf,, k = 
0, 1/2, 1, 3/2, 2, 5/2, 3, respectively, are 



0, 0, 0, 
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Summing up the results for Co( s ) an d C( s ); flA.8j) . (IA.11|) . (1A.14[) . one obtains 
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B1/2 = B\ = B 2 = B 3 = 0. 



(A.15) 



Appendix B. The heat kernel coefficients for a cone with the 5-potential 

In the zeta function (I3.23P we again pick out the term with n = 

(B.l) 



C(a) = CoOO+COO, 
(o(s) = -C(s) J Ayy^-Hl + gl^K^y)}, 

~ °° roc j 

C(s) = C{s)Y,(npy- 2s / d** 1 " 2 ' — ln[l 

n=l dZ 



(zas)] ; f = np. 
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In ((s) we use the uniform asymptotic expansion of the modified Bessel functions 

Mi + ^h^hi^^M, t = ^=L=, (b.2) 



i=l 
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With ( 1B.2j) . the integrals over z converge in the strip (5 — j)/2 < Res < 5/2, where j 
is the power of t, and 
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The contributions from C dlv ( s ) to the heat kernel coefficients we calculate in the linear 
in g approximation by making use of ( 14.21) . These results are presented in the fourth 
column of the table IBU 

In the same approximation Co(0) is given by 



M = 9 -C(s)J™dyy^-^-{I (y)K (y)}. 
The integral (1B.4|) converges in the region 0<Res<l/2, because 

±[I (y)K (y)] = -- - (in | + 7 ) y + . . . , y -> 0; 

d 13 1 

-[I (y)K (y)] = +..., „_>«, 



(B.4) 



(B.5) 



To perform the analytic continuation of (1B.4|) to the domain Re s > 1/2 one has to add 
and subtract from the integrand in ( IB. 41) several terms of its small y asymptote 



Co(s) = £C( a ) 



+ 

where the singular part 
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1 -2s w ~ y 3-2s (3-2s) 2 _ 

has simple poles at the points s = 3/2 and s = 1/2. We use (1B.7j) to calculate the 
contribution of Co( s ) to the heat kernel coefficients B , Bu2, and 5i. These contributions 
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Table Bl. The contributions to the heat kernel coefficients Bk coming from Co(s) and 
The results are obtained in the linear in g approximation. 
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are generated by the infrared singularities of the spectral density in the problem under 
consideration. 

Similarly, the y — > oo asymptote (IB. 51) should be added and subtracted to continue 
the integral ( IB. 41) to the region Re s < 0. In this way we obtain 



f oo 
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(B.8) 



(B.9) 



y 20Fc 2s r(s)r(3/2-s) 
According to (|4.2|) . the singular part (IB.7P gives contribution to the heat kernel 
coefficients Bo, By 2, and B\ while singular part ( 1B.9|) contributes to the heat kernel 
coefficients -B3/2, B 2 , and 5 5 / 2 - These are given in the second column of table IBll 

Summing contributions from £ and (, one obtains the heat kernel coefficients for 
a scalar massless field considered on a cone with semitransparent boundary conditions. 
These results are presented in the last column of table | 



References 

[1] Kibble T W B 1976 J. Phys. A: Math. Gen. 9 1387 

Kibble T W B 1980 Phys. Rep. 67 183 
[2] Vilenkin A 1985 Phys. Rep. 121 263 

Vilenkin A and Shellard EPS 1994 Cosmic Strings and other Topological Defects (Cambridge: 
Cambridge University Press) 

Hindmarsh M and Kibble T W B 1995 Rep. Prog. Phys. 58 477 
[3] Bcvis N, Hindmarsh M, Kunz M and Urrestilla J 2008 Phys. Rev. Lett. 100 021301 
[4] Pogosian L, Tye S-H H, Wasserman I and Wyman M 2009 J. Cosmol. Astropart. Phys. 0902 013 
[5] Brandenberger R, Firouzjahi H, Karouby J and Khosravi S 2009 J. Cosmol. Astropart. Phys. 0901 
008 



Vacuum energy in conical space 



29 



[6] Masahiro Kawasaki, Koichi Miyamoto, Kazunori Nakayama, Gravitational waves from kinks on 

infinite cosmic strings. farXiv:1 002.0652 
[7] Berezinski V, Hnatyk B and Vilenkin A 2001 Phys. Rev. D 64 043004 

[8] Brandenberger R, Cai Y F, Xue W and Zhang X M 2009 Cosmic ray positrons from cosmic strings, 
larXiv:0901. 34741 

[9] Chernoff D F and Tye S-H H 2007 Cosmic string detection via microlensing of stars, 
larXi^:0709.1139l 
[10] Sarangi S and Tye S-H H 2002 Phys. Lett. B 536 185 

Copeland E J, Myers R C and Polchinski J 2004 J. High Energy Phys. JHEP06(2004)013 
Dvali G and Vilenkin A 2004 J. Cosmol. Astropart. Phys. 03 010 
[11] Kibble T W B 2004 Cosmic strings reborn? |arXiv:astro-ph/0410073| 

[12] Davis A-C and Kibble T W B 2005 z Fundamental cosmic strings Contemp. Phys. 46 313 

[13] Polchinski J 2005 Cosmic Superstrings Revisited AIP Conf. Proc. 743, 331-340 (2005); Int. J. 

Mod. Phys. A 20, 3413-3415 (2005); [arXTv:hep-th/0410082 ^2. 
[14] Sakellariadou M 2010 Formation & evolution of cosmic superstrings: a short review Fortschr. Phys. 

No 2/3 larXiv: 1001 .17521 

[15] Vilenkin A 2005 Cosmic strings: progress and problems In "Inflating Horizons of Particle 
Astrophysics and Cosmology", ed. by H. Suzuki, J. Yokoyama, Y. Suto and K. Sato (Universal 
Academy Press, Tokyo, 2006) ]arXiv:hep-th/0508135| r2 

[16] Copeland E J and Kibble T W B 2010 Proc. R. Soc. A 466 623 

Copeland E J, Pogosian L and Vachaspati T 2011 Seeking String Theory in the Cosmos, 
larXnF:1105.0207l 

[17] Vilenkin A 1981 Phys. Rev. D 23 852 

[18] Hiscock W A 1985 Phys. Rev. D 31 3288 

[19] Linct B 1985 Gen. Rel. Grav. 17 1109 

[20] Sokolov D D and Starobinsky A A 1977 Dokl. Akad. Nauk SSSR 234 1043-1046 

Sokolov D D and Starobinsky A A 1977 Sov. Phys., Dokl. 22 312-313 (Engl. Transl.) 
[21] Nesterenko V V, Lambiase G and Scarpetta G 2002 Ann. Phys., NY 298 403 
[22] Gott J R 1985 Astrophys. J. 288 422 

[23] Khusnutdinov N R and Bordag M 1999 Phys. Rev. D 59 064017 

[24] Bordag M, Klimchitskaya G L, Mohideen U and Mostepanenko V M 2009 Advances in the Casimir 

Effect (Oxford: Oxford University Press) 
[25] Milton K A 2010 Local and Global Casimir Energies: Divergences, Renormalization, and the 

Coupling to Gravity, invited review paper to Lecture Notes in Physics volume in Casimir physics 

edited by Diego Dalvit, Peter Milonni, David Roberts, and Felipe da Rosa arXiv:1005.0031vl 

[hep-th] 

[26] Nesterenko V V, Pirozhenko I G and Dittrich J 2003 Class. Quantum Grav. 20 431 

[27] Nesterenko V V, Lambiase G and Scarpetta G 2001 J. Math. Phys. 42 1974 

[28] Nesterenko V V 2008 J. Phys. A: Math. Gen. 41 164005 

[29] Whittaker E T 1904 Proc. London Math. Soc. 1 367 

[30] Nisbct A 1955 Proc. R. Soc. London Ser. A 231 250 

[31] Stratton J A 1941 Electromagnetic Theory (New York: McGraw-Hill) 

[32] Nesterenko V V 2006 J. Phys. A: Math. Gen. 39 6609 

[33] Milton K A, Nesterenko A V and Nesterenko V V 1999 Phys. Rev. D 59 105009 
[34] Lambiase G, Nesterenko V V and Bordag M 1999 J. Math. Phys. 40 6254 

[35] Albeverio S, Gesztesy F, H0egh-Krohn R and Holden M 2005 Solvable Models in Quantum 
Mechanics 2nd edn with an Appendix by Pavel Exner, AMS Chelsea Publishing, v 350 (American 
Mathematical Society, Providence, RI, USA) 
[36] Abramowitz M and Stegun I (eds) 1972 Handbook of Mathematical Functions (New York: Dover) 
[37] Bordag M Pirozhenko I G and Nesterenko V V 2005 J. Phys. A: Math. Gen. 38 11027 
[38] Newton R G 2002 Scattering Theory of Waves and Particles (New York: Dover) 



Vacuum energy in conical space 



30 



[39] Nesterenko V V, Lambiase G and Scarpctta G 2004 Calculation of the Casimir energy at zero and 
finite temperature: Some recent results, Rivista del Nuovo Cimento 27 No 6 1-74 
Nesterenko V V and Pirozhenko I G 2000 J. Math. Phys. 41 4521 
Nesterenko V V and Pirozhenko I G 1999 Phys. Rev. D 60 125007 
Klich I and Romeo A 2000 Phys. Lett. B 476 369 

Gradshteyn I S and Ryzhik I M 2000 Table of Integrals, Series and Products 6th edn (New York: 
Academic) 

Cavero-Pelaez I Milton K and Kirsten K 2007 J. Phys. A: Math. Gen. 40 3607 
Scandurra M 2000 J. Phys. A: Math. Gen. 33 5707 

Scandurra M 2000 Vacuum energy in the background of delta potentials Ph.D. Thesis, Advisor: M. 

Bordag |arXiv:hep-th/0011151^ 2 
Bordag M, Nesterenko V V and Pirozhenko I G 2002 Phys. Rev. D 65 045011 
Midler W 1998 Gommun. Math. Phys. 192 309 

Hurt N E 1983 Geometric Quantization in Action (Boston: Reidel) section 0.10 
Fursaev D V 1994 Phys. Lett. B 334 53 
Fursaev D V 1994 Class. Quantum Grav. 11 1431 
Fursaev D V 1997 Class. Quantum Grav. 14 1059 
Brevik I, Ellingsen S A and Milton K A 2009 Phys. Rev. E 79 041120 
Ellingsen S A, Brevik I and Milton K A 2009 Phys. Rev. E 80 021125 
Brevik I, Ellingsen S A and Milton K A 2010 Int. J. Mod. Phys. A 25 2270 
Milton K A, Wagner J and Kirsten K 2009 Phys. Rev. D 80 125028 
Ellingsen S A, Brevik I and Milton K A 2010 Phys. Rev. D 81 065031 
Rezaeian A H and Saharian A A 2002 Class. Quantum Grav. 19 3625 
Bezerra de Mello E R Bezerra V B and Saharian A A 2007 Phys. Lett. B 645 245 
Saharian A A 2007 Eur. Phys. J. C 52 721 

Bezerra de Mello E R, Bezerra V B, Saharian A A and Tarloyan A S 2006 Phys. Rev. D 74 025017 
Bezerra de Mello E R and Saharian A A 2009 J. High Energy Phys. JHEP04(2009)046 
Bordag M and Pirozhenko I G 2001 Phys. Rev. D 64 025019 



